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ABSTRACT 

We  find  periodic  /-band  variability  in  two  ultracool  dwarfs,  TVLM  513-46546  and  2MASS  J00361617  +  1821104, 
on  either  side  of  the  M/L  dwarf  boundary.  Both  of  these  targets  are  short-period  radio  transients,  with  the  detected 
/-band  periods  matching  those  found  at  radio  wavelengths  (P  =  1.96  hr  for  TVLM  513-46546  and  P  =  3  hr  for 
2MASS  J00361617  + 1821104).  We  attribute  the  detected  /-band  periodicities  to  the  periods  of  rotation  of  the  dwarfs, 
supported  by  radius  estimates  and  measured  v  sin  i  values  for  the  objects.  Based  on  the  detected  period  of  rotation 
of  TVLM  513-46546  (M9)  in  the  /  band,  along  with  confirmation  of  strong  magnetic  fields  from  recent  radio 
observations,  we  argue  for  magnetically  induced  spots  as  the  cause  of  this  periodic  variability.  The  /-band  rotational 
modulation  of  the  L3.5  dwarf  2MASS  J00361617  + 1821 104  appeared  to  vary  in  amplitude  with  time.  We  conclude 
that  the  most  likely  cause  of  the  /-band  variability  for  this  object  is  magnetic  spots,  possibly  coupled  with  time- 
evolving  features  such  as  dust  clouds. 

Subject  headings:  methods:  data  analysis  —  stars:  low-mass,  brown  dwarfs  —  stars:  rotation  —  stars:  spots  — 
stars:  variables:  other  —  techniques:  photometric 


1.  INTRODUCTION 

The  term  ultracool  dwarf  encompasses  very  low  mass  stars 
and  brown  dwarfs,  covering  spectral  types  late  M,  L,  and  T. 
As  the  temperature  decreases,  spectral  type  moves  from  M  to 
L,  with  the  characteristic  TiO  and  VO  bands  of  M  dwarf  spectra 
becoming  replaced  by  broad  neutral  alkali  and  iron  hydride 
lines  in  L  dwarfs  at  roughly  2000  K  (Kirkpatrick  et  al.  2000). 
Although  the  change  in  spectral  features  in  progressing  between 
spectral  types  M  and  L  is  well  characterized,  the  processes 
occurring  in  their  cool,  neutral  atmospheres  are  not  fully  un¬ 
derstood.  Also  uncharacterized  are  any  differences  in  the  phys¬ 
ical  nature  or  extent  of  surface  features. 

One  of  the  main  reasons  for  photometric  variability  studies 
of  late  M  and  L  dwarfs  is  to  probe  the  physical  nature  of  their 
atmospheres.  Ultracool  dwarf  variability  is  generally  attributed 
to  two  sources:  the  presence  of  magnetic  spots  or  dust  clouds 
(Martin  et  al.  2001;  Gelino  et  al.  2002,  hereafter  G02;  Rock- 
enfeller  et  al.  2006,  hereafter  R06).  Surface  features  such  as 
magnetic  spots  or  dust  clouds  may  cause  optical  modulation 
as  the  object  rotates  and  in  certain  cases  allow  a  measurement 
of  the  period  of  rotation  of  a  dwarf.  For  either  of  these  scenarios 
to  be  accepted  as  the  dominant  source  of  photometric  variability 
for  a  particular  ultracool  dwarf,  they  must  explain  the  type  of 
variability  observed,  whether  it  is  aperiodic  or  periodic.  It  is 
sometimes  possible  to  constrain  the  probable  cause  of  vari¬ 
ability  for  an  object,  by  comparing  multiwavelength  photom¬ 
etry  to  photometric  signatures  generated  using  atmospheric 
models  such  as  those  of  Allard  et  al.  (2001).  This  method  has 
recently  been  used  to  indicate  magnetically  induced  cool  spots 
as  the  most  likely  cause  of  periodic  variability  ( P  =  3.65  ± 
0.1  hr)  in  the  M9  V  dwarf  2MASS  1707+64  (R06). 
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There  are  indications  of  changes  in  the  physical  nature  or 
the  extent  of  surface  features  when  moving  from  spectral  type 
M  to  L.  Time-resolved  photometry  of  a  statistical  sample  of 
M  and  L  dwarfs  in  directly  comparable  data  sets  suggests  that 
variability  is  more  common  in  L  dwarfs  than  late  M  dwarfs 
(Bailer-Jones  &  Mundt  2001,  hereafter  BJM01;  R06).  Other 
photometric  studies  have  demonstrated  a  high  level  of  L  dwarf 
variability,  with  a  large-scale  L  dwarf  survey  reporting  vari¬ 
ability  in  seven  out  of  18  L  dwarfs  at  the  95.4%  confidence 
level,  with  a  further  five  targets  showing  variability  at  confi¬ 
dence  levels  of  ~80%  (G02). 

Another  characteristic  of  L  dwarf  photometric  variability  is  a 
lack  of  stability  in  both  periodicity  and  amplitude.  Aperiodic 
variability  is  the  dominant  type  of  L  dwarf  variability  found  in 
the  surveys  by  BJM01  and  G02.  Two  of  the  seven  variables  in 
the  latter  survey,  2MASS  0746+20AB  and  2MASS  1300+19, 
displayed  significant  peaks  in  a  CLEAN  periodogram  (Roberts 
et  al.  1987)  but  did  not  display  persistent  periodicity  throughout 
the  data  set.  Koen  (2006)  reported  a  2.4  hr  periodicity  for  2MASS 
0605—22342,  which  persisted  over  3  days  but  decreased  sig¬ 
nificantly  in  amplitude  throughout  the  observations.  This  ap¬ 
parent  paucity  of  stable  periodic  variability  may  be  linked  to 
spectral  type,  although  there  has  been  no  long-timescale  moni¬ 
toring  of  ultracool  dwarf  variability  to  establish  this  statistically. 

It  has  been  suggested  that  rapid  evolution  of  atmospheric 
features  such  as  dust  clouds  is  the  most  likely  cause  of  observed 
aperiodic  variability  in  L  dwarfs  (G02).  However,  it  may  be 
possible  that  magnetic  spots  on  L  dwarfs  also  produce  signif¬ 
icant  periodic  modulation  but  that  evolving  features  such  as 
dust  obscure  the  underlying  period.  One  theory  that  may  ex¬ 
plain  the  characteristics  of  L  dwarf  variability  is  the  masking 
hypothesis,  proposed  by  BJM01,  which  suggests  that  if  the 
timescale  of  evolving  surface  features  is  shorter  than  the  period 
of  rotation  of  an  object,  the  periodic  modulation  of  its  light 
curve  will  be  inhibited.  In  order  to  assess  the  validity  of  such 
theories  of  variability  for  a  particular  object,  constraints  on  the 
value  of  its  period  of  rotation  from  v  sin  i  measurements  and 
radius  estimates  are  extremely  useful.  Limits  on  the  period  of 
rotation  determine  whether  a  period  detected  in  photometric 
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time  series  is  rotational  and  allow  an  evaluation  of  the  stability 
of  the  photometric  periodicity  in  relation  to  it. 

Recent  Very  Large  Array  (VLA)  observations  of  the  M9 
dwarf  TVLM  513-46546  (hereafter  TVLM  513)  provided  evi¬ 
dence  of  a  1.96  hr  stable  periodicity  (Hallinan  et  al.  2006, 
hereafter  H06;  Hallinan  et  al.  2007,  hereafter  H07).  This  pe¬ 
riodicity  was  found  to  be  consistent  with  a  coherent  radiation 
mechanism  generated  via  the  electron  cyclotron  maser  insta¬ 
bility  operating  in  the  low  plasma  density  regions  above  the 
magnetic  poles  of  the  dwarf,  requiring  TVLM  513  to  possess 
extremely  strong  (~kilogauss)  magnetic  fields.  Such  strong 
fields  were  later  confirmed  for  other  ultracool  dwarfs  by  Reiners 
&  Basri  (2007).  H06  also  argues  that  the  L3.5  dwarf  2MASS 
J00361617  + 1821 104  (hereafter  2MASS  J0036+18)  requires 
the  same  mechanism  to  explain  the  properties  of  its  radio  emis¬ 
sion  (Berger  et  al.  2005,  hereafter  B05),  which  are  almost  iden¬ 
tical  to  those  of  TVLM  513. 

The  extremely  strong  magnetic  fields  of  TVLM  513,  and 
most  likely  for  2MASS  J0036  +  18,  might  be  expected  to  pro¬ 
duce  the  necessary  temperature  gradient  to  create  magnetic 
spots,  thus  causing  photometric  modulation.  Consequently,  we 
decided  to  photometrically  monitor  these  dwarfs  on  either  side 
of  the  M/L  boundary  for  /-band  variability. 

2.  DATA  ACQUISITION  AND  OBSERVATIONS 

2.1.  Supporting  Radio  Data 

The  supporting  radio  observations  discussed  in  this  Letter 
involve  ~10  hr  of  VLA  observations  of  TVLM  513  at  a  fre¬ 
quency  of  8.4  GHz  on  2006  May  20  and  ~10  hr  at  a  frequency 
of  4.9  GHz  on  2006  May  21.  These  observations  and  the  ac¬ 
companying  results  are  detailed  in  H07.  The  4.9  GHz  radio 
data  of  2MASS  J0036+18  is  that  of  B05  and  was  obtained 
from  the  VLA  archives  (project  AB1052). 

2.2.  I-Band  Observing  Strategy 

Choices  made  with  regard  to  the  optical  observations  are  as 
follows:  (1)  the  TVLM  513  observations  were  chosen  to  co¬ 
incide  with  the  nights  of  additional  VLA  observations  (H07); 

(2)  due  to  the  very  red  optical  colors  of  these  dwarfs,  the  / 
band  was  chosen  to  monitor  these  objects  for  variability;  and 

(3)  high  time-sampling  was  performed  for  sensitivity  to  the 
rotational  modulation  of  these  rapid  rotators. 

2.3.  TVLM  513  Observations 

The  four-night  photometric  campaign  of  TVLM  513  (in,  ~ 
15.09)  was  carried  out  between  2006  May  18  and  21,  UT.  Ob¬ 
servations  were  performed  on  the  United  States  Naval  Obser¬ 
vatory  (USNO)  1  m  reflector  in  Flagstaff,  Arizona,  from  May 
18  to  20.  A  Tektronix  2048  x  2048  CCD  was  used,  with  a  pixel 
scale  of  0.68"  pixel  ',  giving  a  field  of  view  of  23'  x  23'.  An 
imaging  cadence  of  ~3  minutes  was  achieved  over  the  ~6  hr  of 
observations  on  each  night.  Additional  observations  were  carried 
out  on  May  20  and  21  at  the  1.8  m  Vatican  Advanced  Technology 
Telescope  (VATT)  at  Mount  Graham,  Arizona.  A  back-illumi¬ 
nated  Loral  CCD  was  windowed  to  a  512  x  512  pixel  frame 
(field  of  view  =1.7'  x  1.7')  to  reduce  readout  time.  This  CCD 
showed  evidence  of  fringing  above  the  sky  noise  level  in  the  I 
band,  so  frames  were  dithered  to  allow  removal  of  fringing  din¬ 
ing  the  data  reduction.  An  imaging  cadence  of  ~1  minute  was 
obtained  for  the  May  21  data  set.  Unfortunately,  due  to  telescope 


tracking  problems,  the  data  obtained  on  May  20  were  of  poor 
quality  and  were  not  used  in  the  analysis. 

2.4.  2MASS  J0036+18  Observations 

Photometric  observations  of  2MASS  J0036+18  (m,  ~ 
16.05)  were  carried  out  on  2006  September  19  UT,  using  the 
USNO  1.55  m  reflector  at  Flagstaff,  Arizona.  The  detector  used 
was  the  Tek2K  CCD,  providing  an  image  scale  of  0.33"  pixeU1 
and  an  11'  x  11'  field  of  view.  Over  the  ~7.5  hr  observation 
of  2MASS  J0036  +  18,  well-exposed  images  were  obtained  ev¬ 
ery  ~8  minutes. 

3.  DATA  REDUCTION 

3.1.  Basic  Reduction 

Science  frames  of  TVLM  513  and  2MASS  J0036+18  were 
processed  using  standard  IRAF5/PyRAF  data  reduction  tech¬ 
niques  and  the  USNO  on-site  reduction  pipeline.  Bias  sub¬ 
traction  and  flat-fielding  were  carried  out  on  all  frames.  The 
2MASS  J0036  +  18  data  were  also  linearized  using  the  Tek2K 
linearity  curve.  The  fringing  pattern  of  the  VATT  CCD  was 
found  to  be  stable,  and  its  removal  was  carried  out  by  producing 
a  fringe-correction  frame  from  the  median  of  a  large  number 
of  dithered,  reduced  science  frames. 

3.2.  I-Band  Differential  Photometry 

Differential  photometry  was  carried  out  on  all  data  sets  in 
order  to  achieve  photometric  precision  of  the  order  of  milli- 
magnitudes  through  the  reduction  of  atmospheric  effects.  Ap¬ 
erture  photometry  was  carried  out  on  the  targets  and  a  selection 
of  reference  stars,  chosen  on  the  basis  of  their  stability,  linearity, 
isolation  on  the  frame,  and  having  a  similar  magnitude  to  the 
target.  Apertures  were  chosen  to  provide  the  highest  signal-to- 
noise  ratio  for  the  target.  Standard  differential  photometry  tech¬ 
niques  were  used,  similar  to  those  outlined  in  BJM01,  and  are 
summarized  below.  Time  series  of  relative  magnitudes  (mrel) 
were  calculated  as  follows:  if  Fi  is  the  instrumental  flux  for 
reference  star  i,  with  a  range  of  reference  stars  i  =  1  ...  M. 
then  the  relative  flux  of  the  rth  frame  (Fr)  is 


The  mrel  values  were  calculated  as  follows,  from  correspond¬ 
ing  target  and  relative  magnitudes,  m,  and  mr,  where  F,  is  the 
measured  flux  of  the  target:  mrel  =  m,  —  mr  =  2.5  log  (FJFf). 

3.3.  Error  Estimation 

Assuming  the  reference  stars  were  well  selected,  the  standard 
deviations  of  their  nonvariable  differential  light  curves  may  be 
used  to  make  an  empirical  estimate  of  the  total  photometric 
errors  for  a  target  of  a  certain  magnitude  (BJM01).  This  ap¬ 
proach  may  be  used  to  estimate  both  formal  and  “informal” 
errors,  such  as  those  from  flat-fielding  and  fringing,  which  may 
be  difficult  to  evaluate  accurately.  The  advantage  of  this  strat¬ 
egy  is  that  it  models  the  expected  error  in  a  target  due  to  noise 
alone,  which  may  be  difficult  to  otherwise  isolate  in  a  target 
with  intrinsic  variability.  The  photometric  errors  in  the  indi¬ 
vidual  reference  star  light  curves  were  calculated  by  plotting 

5  Image  Reduction  and  Analysis  Facility. 
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Fig.  1 . — Lomb-Scargle  periodograms  of  raw  radio  and  optical  data.  I- band 
photometry  of  2MASS  J0036+18  and  a  field  star  of  similar  magnitude,  ob¬ 
tained  on  USNO  1.55  m,  2006  September  19,  and  4.9  GHz  total  intensity 
radio  emission  (B05),  are  also  shown. 

their  formal  photometric  errors  from  IRAF  (crIRAF)  against  the 
standard  deviations  of  their  light  curves  (<rrms).  The  free  param¬ 
eters  of  the  first-order  polynomial  fit  to  this  plot  ( a  and  b)  were 
then  used  to  give  the  errors  in  the  ;th  reference  star  light  curve 
(8m):  8m ,  =  a  +  bo IRAF.  The  error  in  the  relative  magnitude 
(8mK])  was  then  calculated  using  the  magnitude  errors  in  the 
target  8m„  and  the  ;th  reference  star,  6m,  (1=1...  M ),  as  dem¬ 
onstrated  by  BJM01: 

(1  ,  2  M 

—)  2  F,2(8mf.  (2) 


3.4.  Period  Determination  and  Significance  of  Variability 

Periodograms  or  power  spectra  were  used  to  search  for  pe¬ 
riodic  variability.  The  Lomb-Scargle  periodogram  (Lomb  1976; 
Scargle  1982)  was  calculated  for  the  /-band  data  of  TVLM  513 
and  2MASS  J0036+18.  The  power  spectra  or  periodograms 
were  searched  for  significant  peaks,  which  signify  periodic  var¬ 
iability.  When  a  significant  peak  was  found,  the  data  were 
phase-folded  to  this  period.  These  phase-folded  light  curves 
were  then  visually  inspected  as  a  final  test — the  scatter  of  the 
light  curves  were  found  to  be  lowest  when  folded  to  the  periods 
detected  in  the  periodograms  and  began  to  show  increased 
scatter  when  folded  at  other  periods.  The  measure  of  variability 
amplitude  used  is  ard,  the  standard  deviation  of  the  target  light 
curve.  However,  the  significance  of  variability  in  a  light  curve 
is  not  simply  the  ratio  orcl/8m,cl:  as  for  a  large  number  of  points, 
statistically  significant  variability  may  be  detected  even  if  erld 
is  only  slightly  larger  than  8mtel  (see  BJM01). 

4.  RESULTS 

Low-amplitude,  quasi-sinusoidal  variability  was  detected  in 
the  /-band  photometry  for  both  objects,  with  periods  identical 
to  those  found  at  radio  wavelengths  (P  =  1 .96  hr  for  TVLM 
513  and  P  =  3  hr  for  2MASS  J0036  +  18);  see  Figure  1  of  this 
Letter  and  Figure  2  of  H07.  These  periods  were  used  with 
radius  estimates  (Dahn  et  al.  2002)  to  calculate  rotational  ve¬ 
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Fig.  2.— -7-band  light  curves  of  2MASS  J0036+ 1 8  and  TVLM  5 13.  (a)  USNO 
1.55  m  raw  2MASS  J0036+18  light  curve,  2006  September  19.  UT,  measured 
using  nine  reference  stars.  ( b )  USNO  1.55  m  raw  light  curve  of  a  field  star  of 
similar  magnitude  to  2MASS  J0036+18.  using  the  same  reference  stars, 
(c)  Binned,  phase-folded  light  curves  of  TVLM  513  data,  obtained  at  the  VATT 
1.8  m,  2006  May  21  (black),  and  the  USNO  1  m,  2006  May  18  (gray),  both 
measured  using  two  reference  stars.  The  error  bars  are  taken  as  Smrel,  the  average 
expected  error  in  the  light  curves. 

locities  of  ~60  km  s_1  for  TVLM  513  and  ~37  km  s_1  for 
2MASS  J0036  +  18.  The  calculated  velocities  are  consistent 
with  v  sin  i  measurements  of  ~60  km  s  for  TVLM  513  (Basri 
2001)  and  36  ±  2.7  km  s  1  for  2MASS  J0036  +  18  (Zapatero 
Osorio  et  al.  2006),  indicating  a  high  inclination  angle,  i.  On 
the  basis  of  these  results,  we  attribute  the  /-band  periods  to  the 
periods  of  rotation,  with  rotation  axes  perpendicular  to  the  line 
of  sight. 

The  detected  rotational  modulation  in  the  TVLM  513  /-band 
data  (arel  =  0.0076,  8mTts ,  =  0.0054),  along  with  confirmation 
of  strong  magnetic  fields  in  the  radio  data,  suggests  that  mag¬ 
netic  spots  are  the  most  probable  cause  of  variability.  In  ad¬ 
dition,  this  periodic  modulation  appeared  to  be  sustained  over 
the  four  nights  of  the  observing  campaign,  being  detected  in 
periodograms  of  the  optical  data  on  all  nights  between  May  18 
and  20  (USNO)  and  in  the  May  21  data  (VATT).  There  is  a 
high  degree  of  correlation  in  the  light  curves  from  both  tele¬ 
scopes,  despite  instrumental  differences  (see  Fig.  2). 

The  rotational  modulation  found  for  2MASS  J0036  +  18  is 
interesting  (arel  =  0.015,  8mre]  =  0.0089),  as  the  amplitude  of 
the  modulation  apparently  changes  within  the  timescale  of  the 
observation,  a  phenomenon  that  is  not  reflected  in  the  light 
curves  of  reference  stars  or  of  a  field  star  of  similar  magnitude 
(see  Fig.  2).  Just  less  than  three  cycles  are  visible  in  the  2MASS 
J0036  +  18  light  curve,  with  the  amplitude  of  the  first  “peak” 
most  closely  matching  the  third  (Fig.  2).  Due  to  the  short  ob¬ 
servation,  the  Lomb-Scargle  periodogram  assumes  that  this  be¬ 
havior  will  repeat  ad  infinitum,  and  thus  a  second  peak  is 
produced  in  the  periodogram  at  a  frequency  of  half  the  fun¬ 
damental  (Fig.  1).  As  the  3  hr  periodicity  provides  a  rotational 
velocity  consistent  with  v  sin  i  measurements,  and  there  is  no 
physical  reason  for  the  6  hr  periodicity,  the  variation  in  the 
peak-to-peak  amplitude  is  most  likely  due  to  some  underlying 
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aperiodic  variability.  If  this  object  is  indeed  an  analog  of  TVLM 
513  as  argued  by  H06,  it  must  have  strong  magnetic  fields. 
Therefore,  we  suggest  that  magnetic  spots  are  the  most  likely 
source  of  the  3  hr  rotational  modulation,  while  rapidly  evolving 
features  may  cause  the  apparent  change  in  amplitude  of  the 
variability.  Such  features  might  include  dust  clouds  or  magnetic 
spots  that  are  changing  in  size  or  temperature. 

5.  DISCUSSION  AND  CONCLUSIONS 

We  have  detected  periodic  modulation  in  the  /-band  pho¬ 
tometry  of  TVLM  513  and  2MASS  J0036+18,  which  we  have 
determined  to  be  rotational  using  radius  estimates  and  v  sin  i 
measurements.  Radio  observations  have  provided  confirmation 
of  strong  magnetic  fields  for  TVLM  513,  and  there  is  evidence 
for  their  existence  on  2MASS  J0036  +  18.  Based  on  the  optical 
and  radio  results,  we  argue  that  magnetic  spots  are  the  cause 
of  the  variability  of  both  these  objects. 

In  the  case  of  2MASS  J0036+18  (L3.5),  this  is  contrary  to 
expectations  that  dust  is  the  cause  of  all  L  dwarf  variability. 
It  has  been  suggested  that  the  fractional  ionization  of  the  pho- 
tospheric  and  atmospheric  plasma  of  ultracool  dwarfs  is  too 
low  to  sustain  magnetic  spot  structures  (G02;  Mohanty  et  al. 
2002).  However,  if  the  magnetic  field  is  large-scale,  as  has  been 
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confirmed  to  be  the  case  for  TVLM  513  (H07),  it  may  produce 
magnetic  spots.  A  plausible  explanation  for  the  2MASS 
J0036+18  result  is  that  the  amplitude  of  the  underlying  rota¬ 
tional  signal  due  to  magnetic  spots  is  itself  modulated  by  evolv¬ 
ing  features  such  as  dust  clouds.  However,  long-timescale  ob¬ 
servations  are  necessary  to  rule  out  random  noise  effects  and 
to  investigate  the  possible  influence  of  evolving  features. 

It  should  be  noted  that  our  observations  were  designed  to 
be  sensitive  to  the  rotational  modulation  of  rapid  rotators,  with 
high  temporal  resolution  photometry  obtained  for  both  objects. 
This  strategy  was  found  to  be  successful,  and  our  next  challenge 
is  to  monitor  2MASS  J0036+18  to  investigate  timescales  for 
possible  evolving  features  by  isolating  their  effects  from  the 
underlying  rotational  modulation. 
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